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Abstract. - We show that functions of type X n = P[Z n ], where P[t] is a periodic func- 
tion and Z is a generic real number, can produce sequences such that any string of values 
Xs, X s +i, Xs+m. is deterministically independent of past and future values. There are no 
correlations between any values of the sequence. We show that this kind of dynamics can 
be generated using a recently constructed optical device composed of several Mach-Zehnder 
interferometers. Quasiperiodic signals can be transformed into random dynamics using nonlin- 
ear circuits. We present the results of real experiments with nonlinear circuits that simulate 
exponential and sine functions. 



Recent experiments with electronic circuits have shown the possibility of communication 
with chaos [1—4]. The interesting question of communication with chaotic lasers has also 
been discussed in [5]. The fast dynamics displayed by optical systems offers the possibility of 
communication at bandwidths of hundreds of megahertz or higher. Very recently, Umeno et 
al. [6] have proposed an optical device implementation of chaotic maps. They rightly claim 
that the development of secure fiber-optic communication systems can have a large impact on 
future telecommunications. One problem in this area is constructing all-optical devices for the 
transmission of high bit rate signals with the appropriate security. Umeno et al. introduce 
multi Mach-Zehnder (MZ) interferometers which implement a very nice class of chaotic maps. 
Other papers have shown that even chaotic communication systems can be cracked if the chaos 
is predictable [7, 8] . 

In the present Letter we will show that using the same experimental setup of Ref. [6] with 
some small modifications and also other physical systems, it is possible to construct random 
maps that generate completely unpredictable dynamics. 

S. Ulam and J. von Neumann [9, 10] proved that the logistic map X n+ i — 4X„(1 — X n ) 
can be solved using the explicit function X n — sin 2 [#7r2™] . Other chaotic maps are solvable 
exactly using, e.g., the functions X n = sin 2 [6>7rfc™], X n = cos[#7rfc n ], and other functions of 
type X n = P[k n ], where k is an integer [11-14]. For instance, X n = sin 2 [#7r3 n ] is the exact 
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Fig. 1 - First-return maps produced by functions X n — sin 2 [#7r2 n ] (a) and X n — sin 2 [#7r3 n ] (b). 



general solution to the cubic map X n+ i = X n (3 — AX n ) 2 . The first-return maps for the 
dynamics of these systems can be observed in Fig. 1. 

In Ref. [6] an optical circuit composed of N MZ interferometers is presented. The scheme 
of this circuit is shown in Fig. 2. In this experimental setup, the input signal is divided into 
N with equal power after passing an 1 x N coupler. The intensity of light is measured at the 
output of the n-th MZ interferometer using a power meter. Each intensity is defined as X n . 
The path length difference is given by AL(n). The values of AL(n) satisfy the relationship 



AL(n + 1) = mAL(n), 



(1) 



where m is an integer, m > 1. 

The transfer function at the output part of the n-th MZ interferometer is given by the 
equation X n — sin 2 [it AL(n)r / X] , where A is the wavelength of the light source and r is the 
effective refractive index of the optical paths of the MZ interferometers. Thus the output 
powers Xi,X2, ■ ■ ■ ,Xn satisfy the equation X n = sin 2 [jTrAL(l)m n /(Xm)]. 

The authors of Ref. [6] proposed to change the initial conditions by changing the wavelength 
of the light source. Umeno et al. [6] performed experiments for the cases m = 2 and m = 3, 
and N — 2. They measured the values X\ and Xi several times changing the wavelength 
in the range from 1560.00 nm to 1560.5 nm. The figures they obtained are approximately 
equivalent to the first-return maps of the logistic and the cubic maps (See Fig. 1.) 

We should stress here that any chaotic map of type X n +\ = f(X n ) is predictable in the 
short term because X n+ i is always defined as a function of the previous value. For instance, 
for the logistic map, whenever X n w 0.1, X n+i 0.36. 




Fig. 2 - Scheme of an experiment described in Ref. [6] with iV Mach-Zehnder interferometers. In this 
experiment a chaotic sequence is generated. 
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In the present Letter, we will show that functions of type 

X n =P{8TZ n ], 



(2) 



where P[t] is a periodic function, T is the period of P[t], 9 is a real parameter and Z is a 
noninteger number, can generate random dynamics in the sense that X n+1 is not determined 
by any string X , X\, . . . , X n of previous values. 

Moreover we will show that the sequence of values X n is such that past values cannot be 
used to predict future values and future values cannot be used to "predict" past values. Fur- 
thermore, for irrational Z, there are no correlations at all between the values of the sequences. 

Let us define the family of sequences 



X 



k.rn.s 



p 



q" l k) [ q - 



(3) 



where k, m and s are integer. The parameters k distinguishes the different sequences. For all 
sequences parametrized by k, the strings of m + 1 values X s , X s+ll X s+2 , . ■ . , X s+m are the 
same. This is so because X„' m ' s = P [T9o(q/p) s (p/q) n ), for all s < n < m + s. So we can have 
an infinite number of sequences that share the same string of m + 1 values. Nevertheless, the 
next value X^+m+i = P [^o(p/<z) m+1 + Tkp m+1 /q\ is uncertain. In general, X s !^ +1 can 
take q different values. In addition, the value X 8 -i, (X s -i = P \T6o(q/p) + Tkq m+1 /p\), is 



s+m . There can be p 
the future 



.X 



s+m i 



also undetermined from the values of the string X s , X s+ i, X s+2 , ■ ■ ■ ,X 
different possible values for X s _i. Thus, for any string X s , X s+1 , X s+2 , ■ 
and the past are both uncertain. In the case of a generic irrational Z, there are infinite 
possibilities for the future and the past. 

We should remark that the functions (J2J, with a noninteger Z = p/q, are not in the 
class X n +i = f(X n ) nor X n +\ = f(X n ,...,X n -m+i), however they can be expressed as 
random maps of type X n +i = f{X n ,I n ), where /„ is a time-dependent random variable (See 
Refs. [15, 16]). That is, the fact that Z is not integer leads to a sort of time dependence. 

Some properties of function X n = sin 2 [OirZ n ] , which is a particular case of , have been 
already studied in previous papers (See e.g. [15-19].) 

Here we will show that there are no statistical correlations between X n and X m (where 
n 7^ to). We will investigate the functions U n = cos(6*7rZ™), which possess zero mean. Note 
that X n = 1 — U%- The values of U n are found in the interval —l<U n <l. Let us define the 
r-order correlations [20,21]: 



E(U nl U n2 ■ ■ ■ U nr ) := / dU [ P (U )U nl U n2 ■■■U n 



(4) 



The invariant density p(U) is given by p(U) — — U 2 ) and Uq = cos(tt9). We have the 

following formula for the correlation functions E(U ni U n2 ■ ■ - U n 



10 

tain 



e id + e -i9\ we ob _ 



dp [cos(6»7rZ ni ) cos(6»7rZ" 2 ) • ■ ■ cos^ttZ"-)]. Considering that cos6> = \ ( 

E(U ni U n2 ■ ■ ■ U nr ) = 2- r J2 5{oiZ ni + o 2 Z'^ + ■■■ + a r Z n ",0), (5) 



where is the summation over all possible configurations (01,02, ... ,cr r ), with a — ±1, 
and 5(n, to) = 1, if n = m or 8{n, to) = 0, if n 7^ to. We will have non-zero correlations only 
for the sets (n 1? n 2l . . . , n r ) that satisfy the equation 



(6) 
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where Cj = ±1. 

It is easy to see that, for Z > 1, 

E(U n U m ) = (7) 

if n ^ m. As a particular case E(U n U n +i) = 0. We also wish to show that the correlation 
functions E(U^U^ +1 ) are zero when i is even and j is odd (or vice versa) i.e. 

E(W n U 2 n l\ 1 )=0. (8) 

So let us calculate E(UjjJ^): E^U 2 ^ 1 ) = ± £ CT S(a 1 Z n + ■■■ + a t Z n + a t+l Z n+1 + 

■ ■ • + a l+2j+1 Z n+1 ,0). Note that E^JJ^) is not zero only if (en + 02 + h f;)Z n + 

(f i+i + o"i + 2 + ■ • • + (Ti + 2j+i)'^™ +1 = 0. For irrational Z, this equation has not solutions. 

Now we wish to consider all the possible correlations E(U ni U n2 ■ ■ ■ U nr ). We should note 
that E can be non-zero in some "trivial" cases (for instance E(U%>), E{UlW^ +1 )) which 
are related to moments E(U^). This does not affect randomness [20,21]. In the language 
of equations X)T=i a iZ rii — 0, this can happen only due to trivial cases as the following 

Z" - z n - z n+1 + z n+1 = 0. 

We will show that all the "nontrivial" correlations are zero for our functions. Suppose 
n r = 2j + 1. Then E(U ni U n2 ■ ■ ■ U nr ) is not zero only if there are solutions for the equations 
Y^i—i < J iZ ni = 0- But these equations can be written in the form 

N + N X Z + ■ ■ ■ + N 2]+1 Z 2l+l =0, (9) 

where iVj are integer, and N 2 j+i =/= 0. For transcendent irrational Z this equation is never 
satisfied. Thus the sequences generated by function X n — sin 2 [6'7rZ"] with a transcendent Z 
are completely uncorrelated. 

Different aspects of the predictability problem are used in references [22-25] as a way to 
characterize complexity and to find distinction between noise and chaos in experimental time 
series. 

In Refs. [22-25] several quantities are introduced in order to determine the true character 
of the time series. All these methods have in common that one has to choose certain length 
scale e and a particular embeding dimension m. The mentioned quantities discussed in these 
articles dsiplay different behaviors as the resolution is varied. According to these different 
behaviors one can distinguish chaotic and stochastic dynamics. 

Using the results of Refs. [15-19] and the present paper, it is possible to prove that functions 
l|2"jl can represent different kinds of dynamics: chaotic time series (with integer Z), random 
maps or unpredictable sequences (with noninteger Z) , and completely uncorrelated sequences 
of independent values (with generic irrational Z.) 

Our functions can be investigated analytically and their complexity can be calculated 
exactly using theoretical considerations [15, 16] . So these functions can be used as very suitable 
models in order to check the predictions of Refs. [22-25]. 

Moreover, we can produce long sequences of values using our models and then, we can study 
them as experimental time series in the framework of the methods presented in Refs. [22-25]. 
In fact, we have investigated the asymptotic behavior of the quantities discussed in Ref. [22- 
25], and our results coincide with those obtained in the mentioned papers. Additionaly, 
we have checked the formula K = X8(X) + h for the complexity of random maps of type 
X n+ \ = f(X n ,I n ), where K is the complexity of the system, A is the Lyapunov exponent of 
the map and h is the complexity of I n , and 9(X) is the Heaviside step function [24]. Details 
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Fig. 3 - First-return maps produced by function X n = sin [8nZ n ] with Z — 1.8 (a) and Z = n (b). 



of the applications of our results in the problem of distinguishing a chaotic system from one 
with intrinsic randomness will be presented in a more extended paper. 

The statistical properties of pseudo-random number generators are discussed in Refs. [26]. 
The authors of these papers have noticed that almost all pseudo-random number generators 
calculate a new pseudo-random number A„+i using a recursive formula that depends on 
the preceeding values X n +i = f{X n , X n —i, . . . , X n - r +i). They have found that the failure of 
these generators in different simulations can be attributed to the low entropy of the production 
rule /() conditioned on the statistics of the input values X n , X n ^\ 1 . . . , A„_ r+ i. Besides, 
all these generators have very strong correlations even at the macrostate level used in the 
simulations [26]. 

We agree with these researchers that this approach, based on the properties of the generator 
rule, is more profound than the empirical tests. 

In this same spirit, we should say that the rule (J2J produces a dynamics where the future 
values are not determined by the past values. In fact, they can be completely unconrrelated. 

Now suppose that we have the same experimental setup of Ref. [6] , which is represented 
in Fig. 2 schematically, but the equation for the path length differences AL(n) will satisfy a 
relationship very similar to Eq. Q), that is AL(n + 1) = ZAL(n), with the change that Z is 
not an integer. In this case the sequence X n of measured light intensities will be unpredictable. 

Fig. shows different examples of the dynamics that can be produced by function X n = 
sin 2 [6*7rZ"] with different noninteger Z. 

Using the properties of function (0) and further investigation we can obtain the following 
results. Function X n = P[4>(n)], where P[t] is a periodic function and 4>(n) is a non-periodic 
oscillating function with intermittent intervals of truncated exponential behavior, would pro- 
duce also unpredictable dynamics. 

Furthermore, we can construct functions of type X n — h[<fi(n)], again with very complex 
behavior, where h(t) is a non-invertible function and 4>(n) is, as before, a non-periodic oscil- 
lating function with intermittent intervals of truncated exponential behavior. Some chaotic 
systems can produce the kind of behavior needed for </>(n) . This physical system can be con- 
structed, for example, with circuits: a chaotic circuit and a circuit with a non-invertible I—V 
characteristic [27]. An experiment with this scheme is reported in [28]. 

However the most interesting fact is that we can construct function <p(n) without using 
previously produced chaotic signals. In this case we plan to use as input to the non-linear 
system only regular signals. 

A consequence of our theory is that a time-series constructed using three periodic signals 
can be transformed into an unpredictable dynamics. The theoretical result is that the following 



6 



EUROPHYSICS LETTERS 



quasiperiodic 
. signal 



nonlinear system 
equivalent (o 
transformation ■ 
Y = expM 




Fig. 4 - Experimental setup to produce unpredictable dynamics using a quasiperiodic signal. 



function can be unpredictable: 

X n = P[Aexp[Q{n)]], (10) 

where P[y] is a periodic function (in some cases it can be just a non-invertible function,) and 
Q(n) is a quasiperiodic function represented by the sum of several periodic functions. 
As an illustrating example let us study the following function 

X n = sin[0(n)], (11) 

where 4>(n) = Aexp[Q(n)], Q(n) = Pi(n) + Pi(n) + P 3 {n), P^n) = a(n - kTj, when kT t < 
n < (k + 1)7V Here T 2 /Ti, T 3 /T 2 and T 3 /Ti are irrational numbers. 

Note that the functions Pi(n) are piece-wise linear. Function Q(n) is also piece-wise linear, 
but it is not periodic. On the other hand, function <j)(n) will behave as a non-periodic oscil- 
lating function with intermittent intervals of finite exponential behavior. At these intervals, 
function X n behaves as function X n = sin 2 [9nZ n ] . 

Similar properties can be found in function l|ll(l if Q(n) — a\ sin(u;in) + a 2 sin(o;2^) + 
a3sin(w3n). In fact, functions sin^n) behaves approximately as increasing linear functions 
whenever uj n m 2nk where k is an integer. We have performed real experiments using the 
setup represented in Fig. 0] 

In our experiments, a quasiperiodic time-series was used as input to an electronic circuit 
that simulates an exponential function [29]. The output of the exponential system is taken 
as the input to a non-linear system that simulates the sine-function [30]. Fig. [S] shows an 
example of the dynamics produced by the experiment. Details of the experiment will be 
presented elsewhere in a more extended paper. 

In conclusion, we have shown that functions of type X n = P[0TZ n ], where P[t] is a T— 
periodic function, 9 and Z are real numbers, can generate random dynamics in the sense that 
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Fig. 5 - Dynamics produced by the experiment. 
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any string of values X s ,X s+ ±, . . . ,X s+m is deterministically independent of past and future 
values. Furthermore, there are no correlations whatsoever between the values of the sequence. 

The experimental setup schematically represented in Fig. 2 (See Ref. [6]), where the path- 
length differences in the Mach-Zehnder interferometers satisfy the equation AL(n +1) = 
ZAL(n) (with noninteger Z), can be used to produce unpredictable dynamics. 

We have performed real experiments with systems that are equivalent to the scheme rep- 
resented in Fig. These experiments corroborate our prediction that, using just static 
non-linear systems, a quasiperiodic signal can be transformed into a random signal. 
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